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ABSTRACT 

In this work we incorporate the newest ISO results on the mid-infrared 
spectral-energy-distributions (MIR SEDs) of galaxies into models for the 
number counts and redshift distributions of MIR surveys. A three- component 
model, with empirically determined MIR SED templates of (1) a cirrus/PDR 
component (2) a starburst component and (3) an AGN component, is developed 
for infrared (3-120 /mi) SEDs of galaxies. The model includes a complete IRAS 
25 /mi selected sample of 1406 local galaxies (z < 0.1; Shupe et al. 1998a). 
Results based on these 1406 spectra show that the MIR emission features cause 
significant effects on the redshift dependence of the K-corrections for fluxes in 
the WIRE 25 /im band and ISOCAM 15 /im band. This in turn will affect deep 
counts and redshift distributions in these two bands, as shown by the predictions 
of two evolution models (a luminosity evolution model with L oc (1 + z) 3 and a 
density evolution model with p oc (1 + z) 4 ). The dips-and-bumps on curves of 
MIR number counts, caused by the emission features, should be useful indicators 
of evolution mode. The strong emission features at ~ 6-8 /mi will help the 
detections of relatively high redshift (z ~ 2) galaxies in MIR surveys. On the 
other hand, determinations of the evolutionary rate based on the slope of source 
counts, and studies on the large scale structures using the redshift distribution 
of MIR sources, will have to treat the effects of the MIR emission features 
carefully. We have also estimated a 15 /im local luminosity function from the 
predicted 15 /im fluxes of the 1406 galaxies using the bivariate (15 /mi vs. 25 /im 
luminosities) method. This luminosity function will improve our understanding 
of the ISOCAM 15 /mi surveys. 
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galaxies 

1. Introduction 

Current research on cosmology is focused on the debate whether the star formation 
history of the universe has a sharp peak in the redshift range of z ~ 1 - 2 (Madau et 
al. 1996), as derived from the Hubble Space Telescope deep surveys (Williams et al. 
1996; Griffiths et al. 1994) and from ground based redshift surveys (Lilly et al. 1995; 
Koo and Kron 1992), or the universe had been forming stars in a relatively constant rate 
since an early galaxy formation time (z ~ 5) until z ~ 1 (Franceschini et al. 1997; 
Rowan- Robinson et al. 1997), as suggested by the cosmic far-infrared (FIR) background 
studies (Puget et al. 1996; Guiderdoni et al. 1997) and studies on the metal abundance in 
Intra-Cluster Medium (Mushotzky and Loewenstein 1997). 

The problem with the HST surveys and ground based redshift surveys is the dust 
extinction, which may hide much of the star formation in early universe from these 
UV/optical surveys. Dust extinction has little effect on mid-infrared (MIR, 4 /iin — 40 jum) 
surveys. Compared to the FIR (>40 fj,m) deep surveys, the MIR surveys also have the 
advantage of using larger detector arrays and having better resolutions; the latter is of 
particular importance because in these bands the depth of a survey is ultimately constrained 
by the confusion limit. Therefore, MIR surveys will play important roles in galaxy evolution 
studies. Indeed, currently several deep MIR surveys are either being carried out using ISO 
(Franceschini et al. 1996; Oliver et al. 1997), or soon to be launched such as WIRE 
(Hacking et al. 1996; Schember et al. 1996) and SIRTF (Cruikshank and Werner 1997). 

On the other hand, it was known before ISO that in the wavelength range between 
3-20 /iin there are several broad band features (see Puget and Leger 1989 for a review), 
which may contribute substantially to the IRAS 12 /im flux of galaxies (Gillett et al. 1975; 
Phillips et al. 1984; Giard et al. 1989; Helou et al. 1991). The best candidates for 
the emitters of these features are some large molecules (~1 nm) of Polycyclic Aromatic 
Hydrocarbon (Leger and Puget 1984; Allamandola et al. 1985), and the features are often 
called PAH features. New ISO observations (Lu et al. 1996, 1998; Vigroux et al. 1996; 
Boulade et al. 1996) show further that these features are present in the MIR spectra of 
most galaxies, and their equivalent widths can be as large as ~10 /im, comparable to or 
even larger than typical band passes for MIR filters (e.g. the band pass of the ISOCAM 
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15 /mi filter is ~6 jim). This indicates that when any of these features moves in or out 
of the band pass of a MIR survey, because of the redshifts of galaxies, very significant 
K-corrections will occur. Interpretations of MIR survey counts have to take this fact into 
account. A pre-ISO attempt on this issue, based on the laboratory PAH spectra (instead of 
the spectra of galaxies), can be found in Maffei (1994). 

It is the primary aim of this paper to develop a new, emiprical model for the MIR 
SED's of galaxies, and incorprate it into the model for number counts of MIR surveys. 
The SED model is built upon available MIR spectra of galaxies taken from published ISO 
observations, and is applied to the 1406 galaxies (with z < 0.1) of the 25 /iin sample of 
Shupe et al. (1998a), selected from the IRAS Faint Point Source Catalog (Moshir et al. 
1992). For each galaxy in the sample a model spectrum is obtained, which has the MIR 
emission features and which can reproduce the four observed IRAS fluxes. These 1406 
spectra are used to determine the K-corrections and their dispersions in our model for 
predictions on number counts and redshift distributions of MIR surveys. The effects of the 
emission features on these predictions are investigated. Detailed comparisons of the model 
results with some ISOCAM MIR surveys (Oliver et al. 1997; Clements et al. 1998) and 
predictions for WIRE survey will be found in a separate paper (Hacking et al. 1998). The 
computer code of the model presented here is available upon request to CX. 

The paper is arranged as follows: the SED model is developed in Section 2. A local 
luminosity function at 15 /im, calculated from the predicted 15 /im fluxes of the 1406 
galaxies, is presented in section 3. In Section 4 we present the model for source counts and 
investigate the effects of the emission features on the MIR number counts. Section 5 is 
devoted to the discussion and Section 6 to the conclusion. 



2. An empirical model for IR SEDs of galaxies 
2.1. Properties of MIR SEDs of galaxies 

The following properties of MIR SEDs of galaxies can be found from the literature: 

I. The MIR SEDs of normal galaxies (without AGNs) in the ~10 /im region are not 
smooth. Indeed there has been evidence for this since the observations of Gillett et 
al. (1975) of the 8-13 /im spectra of M 82 and NGC 253, the two nearby prototype 
starburst galaxies. Later ground-based observations (Roche et al. 1991) and the 
IRAS low-resolution spectral observations (Cohen and Klein 1989) show further that 
the 11.3 /xm feature and the 8 /im feature are very common in the 8-13 /im spectra of 



-4- 



'HII' galaxies (i.e. starburst galaxies). Most recently, Lu et al. (1996, 1998) obtained 
MIR spectra (3-11.7 /im) of a sample of late-type galaxies using ISO. These galaxies 
are part of a larger set of star-forming galaxies, selected to cover a broad range of 
infrared and optical properties, in the ISO key project on normal galaxies (Helou et 
al. 1996). The spectra were taken with ISOPHT operated in its PHT-S spectroscopic 
mode (Lemke et al. 1996) with a triangular chopping on both sides of the galaxy 
on the sky. In the wavelength range covered by the observations, 3-11.7 /mi, several 
prominent broad band features (e.g. at 6.2 /im, 7.7 /im, 8.6 /im and 11.3 /mi) are 
clearly visible in almost every spectrum in the sample. ISOCAM CVF spectrum of 
NGC 5195 (the small companion of M51) shows further that beyond the PHT-S long 
wavelength boundary there is still another broad feature at ~12.7 /im (close to the 
[Nil] line at 12.8 /im; Boulade et al. 1996). All of these features are presumably 
emission bands of PAH molecules corresponding to different C-C and C-H modes 
(Leger and Puget 1984; Allamandola et al. 1985). 

2. Viewing through the gallery of ISO MIR spectra of galaxies (Lu et al. 1996, 1998; 
Boulade et al. 1996; Vigroux et al. 1996; Metcalfe et al. 1996; Acosta-Pulido et 
al. 1996), of Galactic star formation regions and reflection nebulae (Verstraete et al. 
1996; Cesarsky et al. 1996a, 1996b; Boulanger et al. 1996), and of Galactic cirrus 
(Mattila et al. 1996), it appears that in general, the PAH features and the underlying 
continuum in the MIR band correlate tightly, indicating that they may have the 
same origin. Except for the HII regions and violent starburst regions in interacting 
galaxies, the shape of the MIR spectra of galaxies, of different Galactic nebulae 
and of Galactic cirrus are remarkably similar. However, in the core of Galactic HII 
regions (Cesarsky et al. 1996b; Verstraete et al. 1996) and in a violent starburst 
region in Arp 244 (the famous Antennae which is a local example of galaxy merger), 
another component of the continuum which increases sharply with wavelength appears 
(Cesarsky et al. 1996b; Vigroux et al. 1996; Verstraete et al. 1996). In the 
framework of the dust model of Desert et al. (1990), this latter component is likely 
to be due to the 3-dimensional very small grains heated by the combined heating of 
multiple photons (Boulanger et al. 1996). The characteristic temperature of this 
component is usually too low for it to be seen at wavelengths <20 /im when the 
heating radiation intensity is below ~ 10 4 times the local ISRF (Cesarsky et al. 
1996b). Noticeably, this component affects mostly the MIR emission at A >12 /im 
and not so much for the emission at shorter wavelengths as shown by the ISO PHT-S 
spectra (3 - 11.7 /im) of galaxies in the sample of Lu et al. (1998) which includes 
both quiescent and starburst galaxies. 

3. From their ground based observations Roche et al. (1991, see also Moorwood 
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1986) found that many active galactic nuclei (AGNs) show distinctively different 
MIR spectra from those of starburst galaxies, in particular the broad band emission 
features are absent. This is confirmed by new ISO observations (e.g. Lutz et al. 
1997a; Genzel et al. 1998). Roche et al. (1991) argued that it is the AGN itself, 
instead of other reasons such as the abundance of certain elements, that affects the 
dust emission process to suppress the MIR features either by alternating the grain 
composition within the nuclear region, perhaps by the selective destruction of grain 
species, or through very different excitation conditions. 



2.2. Our model 

Our model is based on the following assumptions: 

Assumption 1: The MIR SED (3-16 /im) of galaxies has three components, each has 
a fixed SED given by the corresponding template (Figures 1-3): (1) a cirrus/PDR 
(Photo-Dissociation Region) component which is dominant in the MIR SEDs of most of 
normal galaxies, (2) a starburst component due to dust heated by the very intense radiation 
in starburst regions which is likely to be the dominant component for violent starburst 
galaxies such as Arp 220, and (3) an AGN component due to the dust emission associated 
with AGNs. 

The adopted MIR spectrum (3-16 jum) of the cirrus/PDR component is presented in 
Figure 1 (the solid curve). Lu et al. (1998) gives two average spectra, one from a set of 10 
"FIR-cold" galaxies with an IRAS 60 /xm-to-100 fxm flux density ratio, fm/i/ fioo^, between 
0.28 and 0.40; and one from a set of 9 "FIR-warm" galaxies with 0.60 < / 60/i / '/ioo/x < 0-88. 
The relative flux accuracy of these spectra should be good to about 20%. The shape of our 
cirrus/PDR template spectrum at wavelengths shorter than ~12 /im is identical to that of 
the average spectrum for the FIR-cold galaxies in Lu et al. (1998). At longer wavelengths 
a tail is taken from the spectrum of NGC 5195 (Boulade et al. 1996). It can be seen that 
the two spectra agree reasonably well at the wavelengths covered by both. 

In the wavelength range of 8.7-16 /im the spectrum of the starburst component 
(Figure 2) is taken from the ISO-CAM CVF spectrum of the overlapping region (region A) 
in Arp 244 (Vigroux et al. 1996), the famous Antennae. This region is where the two disks 
crash, and where much of the star formation in this merging system is taking place. The 
emission features are still visible in the spectrum, but most of the energy is in the continuum 
which has a relatively steep slope. Shorter than 8.7 /im, the shape of the spectrum of region 
A in Arp 244 is very close to that of the average spectrum for the FIR-warm galaxies in 
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Lu et al. (1998). Therefore shorter than 8.7 /im the starburst component is taken from the 
latter (Fig. 2) which extends to 3 /im. Interestingly, longer than 8.7 /im the MIR continuum 
in the average spectrum of the FIR-warm galaxies of Lu et al. is not as strong as in the 
starburst region in Arp 244, suggesting that in these galaxies a significant fraction of the 
MIR emission is emitted outside the starburst regions (i.e. associated with the cirrus/PDR 
component). 

It should be noted that in the short wavelength channel of ISOPHT-S (<5 /im) the 
mean spectra of Lu et al. (1998) are rather noisy, and show no significant evidence for 
any detection of features (e.g. the 3.3 feature) at A <5 /im. Thus both cirrus/PDR 
template and the starburst template presented above are approximated by flat spectra in 
the wavelength range 3 /im < A <5 /im, determined by the average flux densities (over the 
same wavelength range) of the mean spectrum of the FIR-cold galaxies and of the mean 
spectrum of the FIR-warm galaxies of Lu et al. (1998), respectively. 

The template of the MIR SED of the AGN component (Figure 3) is based on the MIR 
spectrum of the Seyfert 2 nucleus of NGC 1068 (Lutz et al. 1997a). Because we are only 
interested in the broad band SED, the narrow emission lines in this wavelength range are 
not included. 

Assumption 2: For a given galaxy with measured IRAS fluxes of fa^ and / 60/Ll , the 25 /im 
flux of the AGN component can be determined using the following 'color- method': 



This assumption is based on the fact that galaxies with low fm^l f^b^ ratios have a large 
chance of exhibiting Seyfert activity (de Grijp et al. 1985). In an optical spectroscopic 
survey for a FIR color selected sample of 358 extragalactic sources with 1 < fa^j f 2 5n < 3.72 
(i.e. —1.5 < ^25,60 < 0), Keel et al. (1994) found that more than 60 percent of the sources 
are Active galaxies (including Seyferts, QSOs and BL Lacs). In Figure 4 is plotted the 
/25/i//i2/i vs. fan / fan color-color diagram of the IRAS 25 /im selected sample of Shupe et 
al. (1998a). Sources with known AGNs (identified using NED) are plotted with different 
symbols. Indeed we find that for log(/ 25 ^/ 'fad > -0.4 (i.e. fan/ 'f 2 5n < 2.5) the vast 
majority of the sources are Active galaxies. In particular, most of the QSOs are on the 
right side of this boundary. On the other hand, most of the normal galaxies (i.e. without 
AGN) have log(/ 25/i / 'fad < -0.3 (i.e. fad fan > 5). Interestingly, many galaxies with 
AGNs, in particular most of the LINERs, also have log(fan/ 'fan) < —0.3. It is likely that 
for these sources the IR fluxes are mainly from the host galaxies and the AGN component 
is not significant. 



•AGN _ 
25^ ~~ 




5 — /s0m//25m 

5-2.5 



(f60/i/f25/i < 2.5); 
(2.5 < f 6 o M /f25 M < 5); 

(f 60/^/^25/* > 5). 



(1) 
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Once is determined, the 12 fim flux of the AGN component is derived as follows: 

rAGN = / 0.6 (f 2 A « N < f 25 , and 0.6 f 2 A 5 f < f^) 

12M \/i2 M (otherwise) 1 J 

The standard /i2/x//25/x ratio of 0.6 for the AGN component is taken from that of the Seyfert 
2 nucleus of NGC 1068 (Lutz 1997a). 

Based on the IRAS data, Giuricin et al. (1995) found that Seyfert 2 galaxies tend to 
have steeper continuum spectra between 12 and 25 //m than those of Seyfert 1 galaxies. 
However, given the large beams of IRAS detectors, it is not clear whether this reflects the 
intrinsic difference between the two types of AGN or whether it is due to the difference in 
the star formation activity around Seyfert Is and Seyfert 2s (Maiolino et al. 1995). In this 
work, we choose to use a single template (Figure 3) to model the MIR SED of both types 
of AGN. 



Assumption 3: Define 
and 



f G _ f rAGN (o\ 



f G _ c rAGN (a\ 

The 12 jum fluxes of the starburst component and of the cirrus/PDR component can be 
determined by the color f^lfvi^ through the following formulae: 

flip if 2%^/ fi2^ — 6); 



rSB 
l 12im 



fG itG _ 
6-2 



(2 < fsV 1 ^ < 6 ); (5) 

(fzSfi/ fl2n — 2) 



and 



f PDR _ f G _ f SB (6) 



In this assumption, we assign a significantly higher standard /25 M //i2^ color ratio 
to the starburst component (/^//rlu = 6) than that of the cirrus/PDR component 
(f25^ R /fu^ R = 2), to be consistent with the fast-rising MIR SED of the starburst 
template (Figure 2). This is also consistent with the anti-correlation between fn^/ f^p an d 
fmn/ fwon °f galaxies (Helou 1986) given that the feon/ fwon ratio is a good star formation 
strength indicator (Bothun et al. 1989; Xu and De Zotti 1989). For example Arp 299 
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(NGC 3690/IC 694), a famous local example of galaxy mergers, has f^/ fn^ = 5.8 and 
/eo/i/ /ioom = 0.95. 

Given Equations p], H H, and |, we can calculate the 12 /im fluxes of the three 
components for a given galaxy with measured /i2 M , and /eo^. These 12 /mi fluxes 
provide the absolute calibrations for the MIR SED templates of these components, the sum 
of which is the model MIR SED (3-16 /im) for the galaxy in question. 

Assumption 4 There are no broad band features at wavelengths longer than 16 /im in the 
IR emission of galaxies. 

This assumption is based on the available IR spectra in the literature. In particular the 
ISOSWS spectrum of M 82 (Lutz et al. 1997b), which covers the wavelength range 3-40 /im, 
does not show any broad band features beyond the 12.7 /im feature. Consequently, for a 
given galaxy, the IR radiation spectrum at wavelengths longer than 16 /im is determined in 
our SED model by fitting the 3 IRAS fluxes at 25, 60 and 100 /im (taking into account the 
real responsivity functions of IRAS detectors), and the flux density at 16 /im predicted by 
the sum of the three MIR components, using a simple cubic spline algorithm (obviously 
the continuum so determined joins at 16 /im the model MIR SED). The advantage of this 
empirical model is that it preserves the observed flux measurements of each source and 
allows the inclusion of a large sample of galaxies in a wide range of observed SEDs. 

Figure 5 is a flow chart showing how actually our model proceeds. 

2.3. Tests of the SED model 

In Figure 6 the predicted MIR spectrum of M 82 and the observed ISOSWS spectrum 
of M 82 (Lutz et al. 1997b) are compared. Allowing for the small apertures of ISOSWS 
channels (ranging from 14" x 20" to 17" x 40"), we have multiplied the IRAS fluxes of 
M 82 at 12, 25, 60, and 100 /im by 0.38, 0.5, 0.5 and 0.38, respectively. Admittedly these 
factors are somewhat arbitrary because at these wavelengths no high resolution maps of 
M 82 are available for us to do aperture photometry. So they are chosen to give the best fit. 
Nevertheless they look reasonable. In particular, the two different factors (0.5 for / 2 5^and 
/ 60ai , 0.38 for /i 2jU and /ioo^) can be understood by means of the known anti-correlation in 
the /i2 M //25/i vs. /eo/V/iooAt diagram (Helou 1986), which suggests that /^and / 60/J are 
more concentrated in the starburst nucleus of M 82. The good agreement between the 
predicted and observed MIR continuum is indeed due to the fine tuning of the aperture 
corrections, while the good agreements between the shapes of the MIR features (which are 
independent of the fitting of the continuum) on the predicted and observed spectra are 
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truly remarkable. 

We made a further test to determine whether the model reliably predicts galaxies' 
MIR SEDs. Ashby et al. (1999) observed with ISOPHT-P a sample of 29 FIR bright 
galaxies (f 6 o M > 5.8 Jy) in several MIR wavelength bands, including the PHT-P 16 /im 
band (52" aperture). Our SED model was applied to a subsample of 18 of the Ashby et al. 
galaxies. The subsample was selected based on four criteria: (1) the signal-to-noise ratio of 
the PHT-P 16 /im flux > 3; (2) detection of the galaxy by IRAS in both the 12 and 25 jum 
bands; (3) optical extent < 2', or when in a galaxy pair the combined size of the pair is 
< 2', (4) the distance to the nearest neighbor (of similar optical brightness) > 4'. The last 
two criteria ensure that the PHT-P observations include most of the emission from the 
region enclosed in the IRAS apertures (~ 4'). Convolving the predicted spectra with the 
PHT-P 16 /im responsivity function, the 16 /im fluxes are predicted and compared with the 
observed /i 6m (Figure 7). Given the uncertainty of the observed /i 6/1 (in most case ~ 30%, 
Ashby et al. 1999), the agreement between the predicted and observed fw^s is quite good. 
The model spectra of the same galaxies are plotted in Figure 8. 

2.4. The Sample 

The IRAS 25 /im selected sample (complete down to / 2 5^ = 0.25 Jy) of Shupe et al. 
(1998a) contains 1456 galaxies, all have measured redshifts and measured IRAS 25 /im 
fluxes. In order to have a local sample free from significant evolutionary effects, our sample 
includes only the 1406 galaxies in Shupe et al's sample which have redshifts < 0.1. Most of 
them are also detected by IRAS at 12, 60 and 100 /im. At 12 /im, 330 (out of 1406) galaxies 
are not detected by IRAS and have therefore only upper limits. For and /ioo^, 18 and 
94 galaxies have upper limits, respectively. We assign a flux density f u = 0.8 x upper limit 
for these cases. 

The SED model thus produces for each galaxy in the sample a SED from 3-120 /im. 
As examples, the SEDs of 50 galaxies with 9 < log(z/L^(25/im)/L ) < 10 are plotted in 
Figure 9. 

3. Local luminosity function at 15 /im 

Several MIR surveys (Oliver et al. 1997; Franceschini et al. 1996) are being carried 
out by ISOCAM in the 15 /im band, the most sensitive one for starburst galaxies among 
ISOCAM filters. However, there is no 15 /im local luminosity function (LLF) available in 
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the literature, because there have been no complete samples of galaxies in 15 /tm band. 
Oliver et al. (1997) used the 12 //m LLF of Rush et al. (1993) in the interpretation of 
their 15 /im deep survey. This may be problematic because different galaxies have different 
15-to-12 /im flux ratios. 

We construct a 15 /mi LLF in this section exploiting the model SEDs of the 1406 
galaxies in our sample. First, for each galaxy in the sample, an in-band flux at 15 /im is 
derived by convolving its SED with the responsivity functions of ISOCAM 15 /mi filter 
(13-18 /tm, ISOCAM User's Manual). The monochromatic flux at 15 /mi is derived from 
this by assuming vf v =constant within the bandpass. When there is only an upper limit 
at 12 /tm, instead of multiplying a by factor of 0.8 to get an approximation of the real flux 
as in the previous section, here we input the 12 /tm upper limit into the SED model and 
obtain an upper limit for the 15 /im flux. The 15 /tm LLF can then be calculated from the 
predicted 15 /tm luminosities, and the corresponding upper limits, of these galaxies. 

Since our sample is 25 /tm-flux limited, not 15 /mi-flux limited, the bivariate method is 
used in the calculation of the 15 /im luminosity function. The basic formula is 

W(L 15li > Zi 5/ J = J 6(Li 5jU > h^\L 25lJt ) p(L 25i m) dL 25fJi (7) 

here ^ is the 15 /tm integral luminosity function, p the 25 /tm differential luminosity 
function (Shupe et al. 1998a), and O the conditional probability function for a galaxy 
with given 25 /im luminosity L 25m to have the 15/t luminosity Li 5fl brighter than Z 15jU . The 
integration is over the entire range of the 25 /tm luminosity. The derivative of \I/ gives the 
differential luminosity function: 

f (90 

0(£l5/i) = / ^7 P(^25 M ) d L 25fM ■ (8) 

In practice, the luminosity functions are calculated in the domain of the logarithms of 
the luminosities. Hence the L^'s and the L 25fl 's in above formulae should be replaced by 
the corresponding logarithms. The sample is binned for both log(L 25M ) and log(L 15At ), and 
Eq(|) can be approximated by 

#og(Li 5 „) = ^m'P[M W = Lj] Aj/5i (9) 

j 

where Lj and Aj are the center and the width of the j-th bin of log(L 25/ J, L$ and Si are the 
center and the width of the i-th bin of log(Li5 M ), and is the probability of galaxies in 
the j-th bin of log(L 2 5^) (i.e. Lj — 0.5Aj < /og(L 25At ) < Lj + 0.5Aj) to have log(Li 5/ J in the 
bin Li (Li - 0.5^ < log(L 15M ) <L i J r 0.5^). 
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In order to take into account the information contained in the upper limits, the 
Kaplan-Meier estimator (Kaplan and Meier 1958) is applied in the calculation of Pij and 
its variance. In what follows we give a brief description of the algorithm. 

For the subsample of galaxies in the j-th bin of log(L 2 5 M ), define F;_i as the probability 
function 

F i _ 1 = P[log(L 15/x )<L i + 0.5(J i ] (10) 



and 

So, 

and 



F t = P[log(L 15/ J < U - 0.55,] . (11) 
Pij = Fi-i - Fi (12) 



Var(Pij) = VariFi-i) + Var(Fi) - 2 x Covar{Fi_ u F<) . (13) 

The algorithm for calculating the KM estimator of Fi, Var(Fi-i) and Var(Fi) can be 
found in Feigelson and Nelson (1985), while the KM estimator for the covariance between 
(for Fi > 0) is (Akritas and LaValley 1997) 

Covar{F l _ l ,F t ) = ^Var{F t ) . (14) 

Finally, the variance of is 

Var[<j>(log(L 15fl ) = LJ] = ^{ Var(P itj )p[log(L 25 ,) = L 3 ] + 

j 

PijVar[p(log(L 2 5n) = L j )]}(A j /5,i) 2 . (15) 

The resultant 15 /zm differential luminosity function is given in Table 1 and plotted 
in Figure 10. A least-squares fit to the integral 15 /im luminosity function by a smooth 
function with the form (Yahil et al. 1991) 



LJ V L 

is carried out. The parameters of this fitting function are given in Table 2. The differential 
15 /im luminosity function specified by these parameters is given in Figure 10 by the solid 
line. 
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4. MIR number counts model and effects of emission features 

4.1. Model for number counts 

Our model for galaxy counts as a function of the limiting MIR flux and redshift follows 
the prescription of Condon (1984; see also Hacking et al. 1987; Hacking & Soifer 1991). 
Here we give a brief outline of the formulation. Details can be found in the above references. 

We consider only one set of cosmological parameters, namely H = 75, = 1 (q = 0.5), 
and the cosmological constant A = 0. 

Let p(L) be the local differential luminosity function (number of galaxies in unit 
comoving volume and unit log(L) interval). When the evolutionary effects are taken into 
account, the LF at redshift z is 

p\L,z)=G(z) P y^y, (17) 

where G(z) is the density evolution function and F(z) the luminosity evolution function. 
The monochromatic luminosity in the rest frame is: 

L v = fAvDl x 10 - 4 ^, (18) 

where f v (hereafter f) is the observed flux density, z the redshift, Dl the luminosity 
distance: 

D L = -^(l + z)[l-(l + z)-^} , (19) 
and K the standard K-correction (Lang 1980): 

Jx 2 S(X)R(X) d\ , 
K = 2.5 log 1 + z) + 2.5 log fA2 ^ 1 . / ; i ; \ ' » • (20) 

/ai 2 S{X/{1 + z))R{X) dX 

The derivative of the comoving volume vs. redshift z is 



dV A ^ 2 dD A 

A 

rr n . .m I ■ ( 21 ) 



— AnD\ — 

dz d z 



c 

H {i~zy* 



where Da is the angular distance: 



^ = 2c [1 -(! + .-)-»] 
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By definition the K-correction is a function of the redshift. However, for a given 
bandpass (R(X)), the K correction also depends on the detailed spectrum (S(X)), which 
varies from galaxy to galaxy. Therefore for the model it is more appropriate to define the 
K-correction as a random variable with its probability distribution depending on the redshift 
and on some other observables. Our prescription for the K-correction explicitly includes 
a luminosity dependence because IRAS data have shown that the /i2^//a5/x color (and 
therefore the MIR SED) is a strong function of the MIR luminosity (Fang et al. 1998). In 
this work we choose to approximate the K distribution by a Gaussian: 



-{K-K { 



U(K, L, z) = — exp 2 ^° J (23) 



\2 



2naK 2cr 



where the mean K and the dispersion o~k are dependent on the luminosity L and the 
redshift z. We have also tried a prescription where the K-correction is a function of the 
FIR-to-blue luminosity ratio instead of the MIR luminosity L. No significant changes are 
found in the results for the number counts. 

Given the above definitions, we have: 

T)(L, z, K) d log(L) dz dK = p'(L, z) U(K, L } z)^d log(L) dz dK (24) 

dz 

where r\ is the counts per unit log(L) interval per unit z interval per unit K interval. 
Substituting the variable L by / according to Eq (|i~8D for given z and K, and integrating 
over K, one obtains the prediction for the differential counts per unit log(/) interval and 
per unit z interval: 

/oo rTy 
f/[(L(f,z,K),z] U[K,L(f,z,K),z}— dK . (25) 



4.2. K-corrections 

The narrower the bandpass[|, the more significant the effect of the emission features 
on the K-correction. Here we compare the predicted K-corrections for three bandpasses: 
ISOCAM 12 /iin band, ISOCAM 15 /an band (ISOCAM User's Manual), and WIRE 25 /mi 
band (Shupe et al. 1998b). The responsivity of these filters are plotted in Figure 11. 

For each of the three bands, K-corrections at different z are calculated for each of the 
1406 galaxies in our sample using its model spectrum and Eq(^). Then the sample is 



The width of a bandpass is better defined by logiXi) — log(Xi), which is not affected by the redshift. 
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binned according to the monochromatic luminosity (rest frame) at the effective wavelength 
of the band, and the means and the dispersions of the K-corrections at different z are 
calculated for the subsample of galaxies in each bin, which are used as estimators of the 
K (L, z) and ctk(L, z) in Eq(E3|). 

In Figures 12a, 12b and 12c the mean K as function of redshift z for galaxies with 
different luminosities are plotted for the three filters, respectively, and the dispersions 
are plotted as error bars. And in Figures 12d, 12e and 12f, the dependences of K and 
<jk on luminosity are plotted for z = 0.5, z = 1 and z = 2 for the same filters. The K-z 
curves for WIRE 25 /im band flux demonstrate two dips near z = 1 and z = 2, which are 
due to the MIR emission features at ~ 11-12 /im and at ~ 6-8 /im (see Figures 1 and 2), 
respectively. The peak around z = 1.5 is caused by the trough at ~ 9-10 /im in MIR SEDs 
of galaxies, which may partially be due to the silicate absorption. For the ISOCAM 15 /im 
band flux, the bump caused by this SED trough moves to z ~ 0.5 and the emission features 
at 7-8 /im cause a dip around z — 1. The K-z curves for the ISOCAM 12 /im band (the 
most broad band among the three) are relatively smooth. 

4.3. Number counts 

Through K-corrections, the MIR features affect the number counts and the redshift 
distributions of deep surveys. Figures 13a, 13b and 13c show Euclidean normalized 
differential counts predicted by different models of the cosmological evolution of galaxies 
for ISOCAM 12 /im, ISOCAM 15 /im and WIRE 25 /im bands. The L oc (1 + zf model 
is chosen because it fits the IRAS 60 /im counts (Pearson and Rowan-Robinson 1996) 
well, and the results from the p oc (1 + z) 4 model are presented to highlight the difference 
between a density evolution model and a luminosity evolution model. The latest local 
luminosity functions (LLFs) at 12 /im (Fang et al. 1998) and 25 /im (Shupe et al. 1998a) 
which improve upon the earlier LLFs in these bands (Soifer and Neugebauer 1991; Rush et 
al. 1993), and the 15 /im LLF presented in Section 3 are used in the calculations of model 
predictions. 

In order to demonstrate the effects of the MIR features on the number counts, in 
Figures 13d, 13e and 13f we show the Euclidean normalized differential counts predicted by 
evolution models based on model spectra of galaxies without the MIR features, which are 
calculated using the smooth fits (power-laws) of the templates of the three components in 
Section 2. 

In all of the three bands considered, the effects of the MIR features on the number 
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counts predicted by the density evolution model (p oc (1 + z) A ) are significant. In 
Figures 13d, 13e, and 13f (no MIR features) the curves of predicted number counts by the 
density evolution model are rather flat, while in Figures 13a, 13b, and 13c the corresponding 
curves show dips and bumps caused by the MIR features. For example in Figure 13b 
(15 /im) the curve of number counts predicted by the density evolution model shows 
two shallow bumps which are due to the emission features in the 11-13 /im wavelength 
range and those in the 6-8.5 /im wavelength range, respectively. The prominent bump 
at ~ 0.1 mJy in Figure 13c (25 /im) is due to the features in both the 6-8.5 /im 
wavelength range and the 11-13 fim wavelength range (see the discussion of Figure 14). 
It appears that in all three bands the MIR features cause only slight distortions on the 
number counts predicted by the no evolution models. 

There is a broad peak on the number counts predicted by the luminosity evolution 
model (L oc (1 + z) 3 ) in all six panels including those for models without the emission 
features. Hence this peak is not due to the emission features. Instead it is a combined 
effect caused by other two factors, one related to the evolutionary model and the other to 
the adopted cosmological model (q = 0.5): The strong evolution (L oc (1 + z) 3 ) ensures 
that before the depth of a survey getting close to the 'edge' of the visible universe (i.e. 
z £ 1) the number counts increase much faster than the Euclidean predictions. Then, 
when the flux level becomes so faint that the characteristic redshifts of galaxies are 
significantly larger than unity (for a given flux, the stronger the luminosity evolution, the 
larger the characteristic redshift), the slope of dV/dz (Eq(^)) starts to decrease rapidly 
with increasing z (decreasing flux), and the number counts drop significantly below the 
Euclidean predictions. 

The effects of the MIR features on the number counts predicted by the luminosity 
evolution model look moderate. The peaks in Figure 13a (12 /im) and in Figure 13b 
(15 /im) are slightly narrower than those in Figure 13d and Figure 13e, and the peak in 
Figure 13c (25 /im) is slightly broader than that in Figure 13f. 

In Figures 14a-f, we compare the redshift distributions of sources at the flux level of 
f u = 0.5 mJy, predicted by different evolutionary models. Note that these distributions are 
for galaxies at this flux level, rather than for galaxies in samples brighter than this flux. 
The flux f u = 0.5 mJy is about the location of the peaks on the predicted number counts of 
the luminosity evolution model (Figure 13), and it is close to the flux limits of all ISOCAM 
deep surveys (Rowan-Robinson et al. 1997; Franceschini et al. 1996) and that of the 
planned WIRE survey (Hacking et al. 1996). 

In all the three bands considered, the no evolution model (both with and without 
MIR features) predicts the mode of redshift distribution of galaxies with f v = 0.5 mJy at 
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z ~ 0.2-0.4. This highlights the reason why the effects of the MIR features are so weak 
on the number counts predicted by this model: because the high redshift galaxies never 
contribute significantly to the number counts, the K-correction (through which the MIR 
features affect number counts) is never significant. 

On the other hand, the MIR features add prominent dips and bumps on the redshift 
distributions predicted by the density evolution model. For example the sharp peak 
at z ~ 1 on the dotted line (density evolution model) in Fig. 14c (25/tm) is due to the 
features in in the 11 — 13/im wavelength range. This peak (after Euclidian normalization) 
reaches its maximum at the flux level when the 11 — 13/zm features of the galaxies are 
redshifted into the 25 /mi filter (L* does not change with redshift in the density evolution 
model, hence for a given flux the L± galaxies have the same redshift). This happens at 
/25/im ~ 0.1m Jy, and therefore a strong bump appears at this flux level on the 25 /tin 
number counts predicted by the density evolution model (Fig. 13c). At fainter flux level 
(/25/im ~ 30fj,Jy) the bump is taken over by sources of z ~ 2 for which the strong emission 
features in the 6-8.5 /mi wavelength range are in the bandpass of the WIRE 25/im camera. 

In Figures 14d, 14e and 14f (no MIR features), redshift distributions predicted by the 
luminosity evolution model are much broader than those predicted by the density evolution 
model and by the no evolution model. The reason for this difference is the following: 
When the K-correction is insignificant, the shape of the redshift distribution is essentially 
determined by the product of the luminosity function and the dV/dz (c.f. Eq([25"D). Thus 
the redshift distributions predicted both by the density evolution model and by the no 
evolution model resemble the 'visibility function' of a local sample (e.g. Shupe et al. 
1998a), with a sharp peak centered at the redshift of the galaxies (for the given flux). On 
the other hand, in the luminosity evolution model, the L+ increases rapidly with redshift, 
causing a much broader peak in the redshift distribution (e.g. Figure 14f). It is because 
of this broadness of the redshift distributions that the effects of the MIR features on the 
number counts predicted by the luminosity evolution model (Figure 13) are relatively weak, 
in particular compared to those predicted by the density evolution model (Figure 13): the 
effects are largely smeared out by including galaxies in a wide range of redshifts for a given 
flux. 



5. Discussion 

MIR emission features, if they exist in the IR SEDs of distant galaxies as in those of 
the local galaxies, can cause significant effects on the K-corrections (Figure 12), especially 
when the bandpass is comparable to the widths of the features (e.g. ISOCAM 15 /mi 
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and WIRE 25 fim bands). This affects the predictions on MIR number counts and the 
redshift distributions of different evolution models. The dips-and-bumps caused by the 
MIR features on the curves of the Euclidean normalized number counts will provide very 
useful diagnostic indicators of the type of evolution that is occurring, which may help to 
distinguish the density evolution model from luminosity evolution model when the redshifts 
are not available. On the other hand, the fluctuations in the slope of the number counts 
caused by the emission features should be taken into account in any attempt of determining 
the evolutionary rate from the slope. 

Dramatic features are caused by the MIR features on the redshift distribution, 
suggesting that certain selection effects in redshift space are expected in MIR surveys. This 
may indeed be beneficial rather than annoying. For example, because of the strong emission 
features at ~ 6-8 /iin, the WIRE 25 /im survey may detect a large number of galaxies of 
z^,3. However, in studies on large scale structures based on deep MIR samples, these 
features must be separated from those caused by large scale inhomogeneities. 



6. Summary and Conclusions 

Our main results are summarized as follows: 

1. An empirical, three-component model for the MIR SEDs of galaxies is developed. 
The templates of (1) a cirrus/PDR component, (2) a starburst component and (3) an 
AGN component are constructed from published ISO observations of MIR spectra of 
galaxies. Tests show that the model can reproduce the MIR spectrum of M 82 and 
the ISOPHT-P 16 p,m fluxes of a sample of galaxies quite well. 

2. The model is applied to a sample of 1406 local galaxies (z < 0.1; Shupe et al. 
1988a). Treating the K-correction as a random variable, its probability distribution 
(approximated by Gaussian) as a function of the redshift z and of the monochromatic 
luminosity L v are calculated for three different bandpasses (WIRE 25 jum, ISOCAM 
15 /im and ISOCAM 12 /im) using these 1406 spectra. 

3. Effects of MIR emission features on the MIR source counts and on the redshift 
distributions are studied for the same MIR bands assuming three different evolution 
models (no-evolution, L oc (1 + z) 3 and p oc (1 + z) 4 ). It is found that predictions by 
luminosity evolution model (L oc (1 + z) 3 ) and density evolution model (p oc (1 + z) 4 ) 
for source counts and redshift distributions in the WIRE 25 /im band and the 
ISOCAM 15 /im band are affected by the MIR emission features significantly, showing 
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dips and bumps which can be identified as being caused by specific MIR emission 
features. 

4. We point out that the dips-and-bumps on curves of MIR number counts will be useful 
indicators of evolution mode. The strong emission features at ~ 6-8 /im will help the 
detections (especially in the WIRE 25 fxm survey) of relatively high redshift (z ~ 2) 
galaxies. On the other hand determinations of the evolutionary rate based on the 
slope of source counts, and studies on the large scale structures using the redshift 
distribution of MIR sources will have to treat the effects of the MIR emission features 
carefully. 

5. A 15 jum local luminosity function is calculated from the predicted 15 /zm fluxes of 
the 1406 galaxies using the bivariate (15 /im vs. 25 /im luminosities) method. This 
luminosity function will improve our understanding of the ISOCAM 15 jum surveys. 
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Table 1. Local differential luminosity function at 15 jim. 



\og(vL v (lh^m)/L & ) 


log(0/(Mpc" 3 mag~ 1 )) 


1 a error 


7.7 


-2.23 


-2.82 


8.1 


-2.88 


-3.34 


8.5 


-2.76 


-3.43 


8.9 


-2.88 


-3.67 


9.3 


-3.13 


-4.06 


9.7 


-3.52 


-4.46 


10.1 


-4.18 


-5.05 


10.5 


-5.06 


-5.80 


10.9 


-5.93 


-6.58 


11.3 


-7.00 


-7.46 


11.7 


-7.91 


-7.91 



Table 2. Fitted Parameters of the 15 /im Luminosity Function 



a f3 log(L*/L ) log(C/Mpc- 3 ) 



0.47 ±0.13 2.20 ±0.13 9.62 ±0.15 



-2.93773 ±0.23 
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Fig. 1. — Solid curve : template of the MIR spectrum of the cirrus/PDR component. 
Crosses with error bars : mean ISOPHT-S spectrum of the ten FIR cold galaxies in the 
sample of Luet al. (1998). Dashed curve : ISOCAM CVF spectrum of NGC 5195 (Bouladeet 
al. 1996). 

Fig. 2.— Solid curve : template of the MIR spectrum of the starburst component. 
Crosses with error bars : mean ISOPHT-S spectrum of of the nine FIR warm galaxies in 
the sample of Luet al. (1997). Dotted curve : ISOCAM CVF spectrum of star-burst region 
A in Arp 244 (Vigrouxei al. 1996). 

Fig. 3.— Template of the MIR spectrum of the AGN component, based on the MIR 
spectrum of the Seyfert 2 nucleus of NGC 1068 (Lutzet al. 1997a). Because we are only 
interested in the broad band SED, the narrow emission lines are not included. 

Fig. 4. — The /25^/ fu^ vs. /25/i/ feon color-color diagram of the IRAS 25 //m selected sample 
of Shupeet al. (1998a). Sources with known AGNs (identified through NED) are plotted 
with different symbols. The rest of the sample are plotted by small dots. The dashed line 
marks the adopted boundary beyond which all the IR emission is considered to be due to 
the AGN component. 

Fig. 5. — Running chart of the SED model. 

Fig. 6. — Comparison between predicted (dotted curve) and observed (solid curve, taken 
from Lutzet al. 1997) MIR spectra of M 82. 

Fig. 7. — Comparisons between predicted and observed (ISOPHT) 16 /im fluxes (fw^prt 
and fiQ^, respectively) for galaxies in the sample of Ashbyet al. 

Fig. 8. — Model spectra of galaxies in the sample of Ashbyet al., overlaid by observed flux 
densities (IRAS: 12 /im, 25 /im, 60 jum and 100 /im; ISOPHT: 16 /im). 

Fig. 9. — Model spectra of individual galaxies in the luminosity range 9 < log(L 2 ^ il / L & ) < 
10. 

Fig. 10. — The 15 /im local differential luminosity function (Table 1). The solid line is the 
smooth fit specified by the parameters given in Table 2. 

Fig. 11. — The bandpasses of the ISOCAM 12 /im and 15 /im filters and of the WIRE 
25 /xm camera. 

Fig. 12. — Panel a. The mean Kq and the dispersion ok (shown by the error bars) of 
K-corrections of the ISOCAM 12 /im flux as functions of redshift z for galaxies in three 
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different luminosity bins. Panel b. Same as Panel a. For the ISOCAM 15 /im band. Panel 
c. Same as Panel a. For the WIRE 25 /iin band. Panel d. Dependences of K and ok on 
the luminosity in the ISOCAM 12 /im band, plotted for three given redshift: z = 0.5, z = 1 
and z — 2. Panel e. Same as Panel d. For the ISOCAM 15 jum band. Panel f. Same as 
Panel d. For the WIRE 25 jum band. 

Fig. 13. — Panel a. Euclidean normalized differential counts in the ISOCAM 12 /im band, 
predicted by three different evolution models. Panel b. Same as Panel a. For the ISOCAM 
15 fim band. Panel c. Same as Panel a. For the WIRE 25 /im band. Panel d. Euclidean 
normalized differential counts in the ISOCAM 12 /im band, predicted by evolution models 
otherwise identical to those in Panel a except that K-corrections are estimated from model 
spectra of galaxies without the MIR features. Panel e. Same as Panel d. For the ISOCAM 
15 /im band. Panel f. Same as Panel d. For the WIRE 25 /im band. 

Fig. 14. — Panel a. Redshift distributions of sources with f v = 0.5 mJy in the ISOCAM 
12 /im band predicted by the three evolution models (note that the distributions are for 
galaxies at this flux level, rather than for galaxies in samples brighter than this flux). Panel 
b. Same as Panel a. For the ISOCAM 15 /im band. Panel c. Same as Panel a. For 
the WIRE 25 /im band. Panel d. Redshift distributions of sources with f v = 0.5 mJy in 
the ISOCAM 12 /im band predicted by the evolution models otherwise identical to those in 
Panel a except that K-corrections are estimated from model spectra of galaxies without the 
MIR features. Panel e. Same as Panel d. For the ISOCAM 15 /im band. Panel f. Same 
as Panel d. For the WIRE 25 /im band. 
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